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ABSTRACT 

We present a new luminous star in M 33 located in the nuclear region. The star shows strong 
Fe II and [Fe II] forest, hydrogen emissions in the spectrum, as well as nebular lines. Ti II and 
Sill lines were detected in absorption, their radial velocity shifted by ^ — 30km/s relative 
to emission lines. The star is variable over seven years with 0^5 variations over a year. We 
studied its spectral energy distribution together with five confirmed Luminous Blue Variables 
and VarA in M33 using homogeneous data and methods. We found the star's bolometric 
luminosity to be logL/L© ^ 6.27, a surface temperature of T~ 16000 K and black body 
temperatures of two dust components of T ~ 900 and 420 K. The new star has properties 
intermediate between the LBVs and VarA (probable cool hypergiant). In the same time it has 
a hot photosphere, LB V-like luminosity and an extensive circumstellar material (strong [Call] 
lines). In these seven luminous variables in M 33 we find the total range of the hot component 
luminosities is 1.0 dex, but that of the dust componets is 2.0 dex. We conclude that the dust 
phenomenon in the luminous variables is temporary and variable, and that dust activity may 
follow strong eruptions. 

Key words: Massive stars - luminous blue variables: general - stars: individual: new Lumi- 
nous Variable - optical spectroscopy: stellar population - galaxies: individual - M 33 



1 INTRODUCTION 

Luminous Blue Variables (LBV) are the most massive stars in on e 
of the final stages of evolution (Humphre ys & DaYidson|[T994h . 
They are very important for understanding of the late stages of 
, evolution, production of WR stars, Supernovae, neutron stars and 
black holes. Their role in evolut ion of galaxies is si gnificant. Stel- 
lar evolution calculations (e.g., iMevnet et al.ll2007l) propose clear 
schemes of massive star evolution, howeve r, the seque nce of stel- 
lar conversions is no t well unders tood yet ( Smith & Ow ocki 2006; 
ISmith & Contill2008l : ISmithll2008h . In order to understand the final 
evolution stages of massive stars, we need to discover and study 
more o bjects at the LBV and LBV-r elated stages. 

IHumphrevs & Davidsonl (Il99 4) reviewed LBV stars and their 
data known up to that time. In the galaxy M33 they selected 
four obvious LBV objects - VarB, VarC, Var2, Va r83, and 
one LBV candidate V532 (GR290, "Romano's star", jRomand 
Spectroscopy has confirmed (Fabrika 2000) the LBV sta- 
tus of V 532. In photometric studies ( Kurtev et al.l l260ll) Its typi- 
cal LBV behaviour has been confirmed a s well. In following pa- 
pers (iFabrika et alboOSl : IViotti et al.ll2006l . l2007h the star has been 
studied in more detail. The star VarA in this galaxy had all the 
properties typical for LBV stars (Hubble & Sandage 1953), it was 
one of the visually brightest stars in M 33 with an intermediate F- 
type spectrum. However, the star shown M-type spectrum (unusual 
for classical LBVs) after its huge brightening in 1950-1953, thus 
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it was classified jHumphrevs et al.lll987l : IHumphrevs et al.l [20061 : 
IViotti et al.ll2006h as a cool hypergiant. 

Searches for LBV ca ndidates in M33 have been done 
using various me t hods dNeese et"aL 
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'Calzetd et al' "1995*; 'Fabrika & Sholukhova' '1995'; "Massey et al.| 
1996; Sholukhova et al. 1997; Fabrika & Sholukhova 1995; 
ISholukhova & Fabrikal l200Ql : ICorral & Herrerol 12003). The main 
technique is to search for Ha-emitting sources coinciding with 
the galaxy stars. Selections of early-type stars were preferable, 
however, k nown LBV objects show various observed colors 
JSterken eT a l. 2008). Some works were ori ginally add resse d 
to a search for SS 433-like objects in M33 dNeese et all Il99ll : 
ICalzetti et all 1 19951 : iFabrika & Sholukhoval ll?95'). SS433 has a 
spectrum very similar to that of a late WN star ( Fabrika 2 004h . 

A new principal step has been done by Massev et al] (l2006h 
in their CCD survey (UBVRI, Ha, [SII], [OIII]) of Local 
Group galaxies. The final catalog in M 33 contains 146,622 stars 
up to 23rd magnit ude with photometrical accuracy of 1-2 %. 
iMassev et al describe a list of Ha emission-line stars in Lo- 

cal Group galaxies. They focused on a search for new LBVs based 
on their spectroscopic similarity to known LBVs. Their number of 
known and suspected LBVs in M33 is 37. Some of their candidates 
coincide with previously described LBV candidates in M33 ( Corrail 
1996; Sholukhova et al. 1997, 1999; Sholukhova & Fabrika 2000; 
[Fabrikall200 0). 

We perfor med photometry of p oint-like Ha emitters in the 
M 33 images of iMassev et al.l (120061) and selected stars with V < 
18.5 mag, B — V < 0.4 having Ha excess. We did independent 



2 Valeev et al 




Figure 1. Nuclear region of M 33, the bright blue objec t is the nucleus. T he 
new star (J2000, 013352.42+303909 or N93351 by Mas sev et all ( l2006h ) is 
indicated. The spectral slit covering both the star and another blue candidate 
is shown. The image was obtained as rgb from Ha, V and [OIII]. 
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Figure 2. Light curve of N 93351 in V-band. The data for 1997- 
1998 were ta ken from Macri et al. ( 2001), the data for 2000-2004 from 
IShporer & M azeh ( 2006). A systematic shift -0.2 mag was applied to the 
second data. Photometrical errors are less than 0.03 mag. 

basic calibrated data" and did photometry of all our targets using 
the task phot in IRAF (http://ssc.spitzer.caltech.edu/irac/dh/). All 
five Spitzer observations were averaged to have one photometrical 
measurement for each target. 



aperture photometry in the UBVRI and Ha images for selected 
stars. Compared to automatic PSF photometry (which is t he only 
possi ble way of dealing with a huge number of stars (Masse vet al.l 
aperture photometry is controlled by eye to avoid confusion 
in crowded (broad-band) and complex (Ha) reg ions. In the selec - 
tion we used the same method as discribed by Fa brika et al.l (Il997h . 
With a di stance modulus of 24.9 mag (Bonanos e t al .12006) our list 
dValeev e t al. 2009) covers bright blue supergiants in the galaxy. 
In our follow up spectroscopy of the list stars we found a new lu- 
minous variable star. Here we report a spectral and photometrical 
study of the star and examine its spectral energy distribution to- 
gether with all the confirmed LBVs and Var A in M 33. 



2 OBSERVATIONS 

For spectroscopy w e use d the SCORPIO focal reducer 
dAfanasiev & MoiseevI l2005h of the Russian 6-m telescope. 
We took spectra of the new star on November 18, 2007 with a 
spectral range and resolution of AA3700 — 7200 A and 12 A, on 
January 7, 2008 (AA3900 - 5700 A and 5 A), and on January 
8, 2008 (AA5700 - 7400 A and 5 A). The slit width was l'\ 
and seeing ^ 1.5^^ The spectra were reduced using standard 
procedures. The accuracy of radial velocity calibration is about 
lOkm/s. 

In Fig.[T]we show the central region of M33, The star is lo- 
cated 34.3^^ (160 pc) from the nucleus to the South-East. The faint 
Ha filaments around the star may be related to it. Our slit position 
was selected to cover the second blue star candidate (Fig.[T])- In the 
HST archival broad band images the star is single and point-like. 

For interstellar absorption estimates, both in this star and the 
other six stars in M33 (Var A, VarB, VarC, Var 2, Var 83 and 
V 532), we used Ha/H/3 emission line ratios (Section|4]) measured 
in surrounding nebulae. We used our own spectra of these targets 
taken with the same tel escope and the MPFS instrument (IFU mode 
jAfanasiev et alJl200lh ) in December, 2004 and October, 2005 and 
with SCORPIO (December, 2004 and October, 2008). The spectral 
data were used for independent estimates of the star temperatures. 

In studying the SEDs of the stars, we used published opti- 
cal photometry (Section© and IR archival data from 2MASS and 
Spitzer. The 2MASS survey of M33 is available in J, H, K bands. 
Spitzer observations of M 33 were made using all four bands of the 
Infrared Array Camera (3.6, 4.5, 5.8, and 8.0 ^im). We used "post- 



3 THE NEW LUMINOUS VARIABLE N 93351 

The new sta r N 93351 has a V m agnitude of 16.16, and B-V=0.28, 
U-B=-0.44 (lMassevetaLll2006h . A light curve of the star is pre- 
sented in Fig|2l The star shows a variability for the past 7 years. It 
shows 0™5 variability over a year. 

In Fig[3]we present spectrum of N 93351 taken with a resolu- 
tion of 5 A. The star has a blue continuum and very bright hydrogen 
emissions with broad wings. The most remarkable feature is Fell 
and [Fe II] emission forest. We detect He I, Ti II and Si II absorp- 
tion lines. He I absorption may be partly filled with nebular emis- 
sion, which is seen in the spectral image. We detect NalDl, 2 and 
Call H, K lines. The last lines may be not only interstellar in origin. 
The spectrum presents strong [OII]A3727 line, Balmer absorp- 
tion edge is practically absent (Figg]). [OIII] AA4959, 5007 lines 
are present, they are faint and crowded with bright Fell emissions. 
There are collisionally excited nebular lines of [S II] AA6717, 6731, 
[Nil] AA6548,6583 (and marginally [NII]A5755) in the spec- 
trum. 

The spectrum is quite typical for LBVs. However we ob- 
serve very strong [Call] AA7291, 7323 emissions, which are not 
present in classical LBVs. These lines a re seen in the cool hy- 
pergi ants Var A in M33 and IRC-hl0420 ('Humphrevs etaLl ll987l : 
Hump hrevs et al . 2006) with extensive circumstellar material. We 
conclude therefore, that N 93351 has intermediate spectrum be- 
tween classical hot LBV stars and the cool hypergiants like Var A. 

The radial velocity of Fe II and nebular lin es is ^ — 120km/s, 
which is in accordance with the HI velocity jReakes & NewtonI 
Il978h in this area of the galaxy. We observe a blue shift of ~ 
— 30km/s in all absorption lines relative to hydrogen emission 
lines. This shift may arise in an expanding photosphere of the star. 
Fell and [Fell] line widths correspond to our spectral resolution 
(5 A, measured in sky lines). For the Ha and Hf3 line profiles we 
performed gaussian analysis for broad and narrow components and 
found the component widths (corrected for the instrumental reso- 
lution) of 800km/s and 180km/s in Ha, 950km/s and llOkm/s in 
HP correspondingly. Radial velocities of the narrow components in 
these two lines are about the same, —106 =b 5 km/s, and velocities 
of the broad components are —125 ± 10 km/s. 

Table[T] presents intensities and equivalent widths of main 
lines, which are not blended. The intensities are given in units 
of intensity of H(3. We present also the line intensities corrected 
for interstellar extinction, where we adopted Ay = 1.0 (Sec- 
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Figure 3. Spectra of N 93351 in blue (top) and red (bottom) regions. Some lines are indicated (upper row), the middle row of vertical bars marks Fell lines, 
the bottom row marks [Fell] lines. CCD fringes are notable at wavelengths longer than ~ 6800 A. 




Table 1. Equivalent widths, relative intensities and dereddened (Ay = 
1.0) relative intensities in units of Hp intensity of main not blended lines. 
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Figure 4. Dereddened spectra of N 93351 (middle) together with stan- 
dard stars O 8 1 (HD 225 160, top) and B 9Ib (LSI V P24, bottom) taken from 
IJacobv etal.llT984) . 



tion|4]). To extimate the extinction we measured Ha and YL/3 line 
fluxes in the surrounding nebula at a distance from 5^^ to 13^^ off 
the star along the slit. We found the fluxes per squared arcsec- 
ond F(Ha) = (1.08 ± 0.04) • lO"^"^ ergW sec and F(H/3) = 
(2.5 it 0.1) • 10~^^ erg/cm^ sec. Corresponding line ratio Ha/H/3 = 
4.0-4.5 gives Ay ^ 1.0 - 1.4. The equivalent width of Na ID 1,2 
absorption ( TabIe[T]) confirms th is value of extinction, a rough lin- 
ear relation jBarbon et al.lll99Ql) between the doublet EW and red- 
dening gives Av = 1.2 zb 0.2. 

Fig.|4] presents spectrum of N 93351 taken with resolution 
12 A corrected for the extinction (Ay = 1.0, E(B-V)=0.32) to- 
gether with standard stars 81 (E(B-V)=0.56) and B9Ib (E(B- 
y)=0. 23), whose dereddened spectra were taken from I Jacob v et al.l 
('1984'). Both the spectral types and reddening values of the standard 
stars are considered as approximate ( Jacoby et al. 1984). Compar- 
ing the spectrum slope and intensity of the B aimer edge with those 
in classified stars we estimate roughly the temperature of N 9335 1 's 
photosphere of 13000 - 16000 K. Absorption Ca II K and Ti II lines 
may be formed in cooler extended atmosphere (the wind). 
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4 SPECTRAL ENERGY DISTRIBUTIONS 

To estimate the temperature and luminosity of this star and other 
six M33 stars (VarA, VarB, VarC, Var2, Var83 and V532) we 
study their SEDs. Our main so urces for the data are optical pho- 
tometry bv MassevetalJ(l2006h (October 2000 - September 2001), 
2MASS (December, 1997), Spitzer (January 2004 - June 2005) and 
our spectroscopy (December 2004 and October 2005). We use these 
data, otherwise describe individual cases. We checked all these 
stars in the literature for their variability so as not to use data re- 
lated to different photometrical states. 

We used our own spectra in order to have independent esti- 
mates of a reddening and a temperature of the stars. In the estimates 
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Table 2. Parameters of the studied stars. Columns present the star name, 
adopted values of extinction, derived bolometric luminosity, temperature of 
the stars (Ts), and temperatures of warm (T^^) and cold dust (Tc) compo- 
nents. All temperatures are given in units of 1000 K, lumimosities in solar 
units. Contributions of the dust components (per cent) to the total luminosi- 
ties are given in brackets. The star (Ls) and dust component (Ld) luminosi- 
ties are given in last columns. 
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of the reddening we used Ha/H/3 flux ratios measured in nebulae 
surrounding the stars. In the Case B of a standard gaseous neb- 
ula this ratio is 2.87 (H7/H/3 = 0.47, = 0.26), and stays 
the same to within 10 % in a wide ran ge of gas temperatures and 
densities ( Osterbrock & Ferlandll2006h . Obvious limitations of this 
method are that the nebula may not be connected with the star, and 
the star may have some additional curcumstellar extinction. In two 
stars VarA and VarC there is no clear morphological association 
with the nebula. 

We used a standard method to fit observed SED components to 
black body (BB) spectra with a reddening (we adopt Rv = 3.1). 
Using the best fit SED solutions we tried to satisfy the indepen- 
dently estimated stellar temperatures (first) and the extinction val- 
ues (second). We fitted the observed SEDs of the stars to BB spec- 
tra, suggesting that one hot spectrum is stellar and one or two cooler 
spectra for dust emission. We kept the temperature of cool compo- 
nents to less than 1500 K. We estimated luminosities integrating the 
extinction corrected SED components. 

In some cases we did not use U data, when this band fluxes 
were obviously weak (because of the B aimer absorption). As a rule 
we did not use R, I (and sometimes J) data, because they can con- 
tain a contribution of free-free emission from stellar winds. Fig.|5] 
presents the observed and intrinsic SEDs of the stars as well as the 
best fits of BB models. Derived parameters of the stars are pre- 
sented in Table 2. Below we describe the individual stars in more 
detail. 

N 93351. The star does not show strong variability. Regarding the 
probable temperature and extinction ranges of the star found in Sec- 
tion 0]), we find the best fit with Ay ^ l.O and T^ 16000 K 
(Fig.|5]). In Table 2 we present a second fit version with Ay = 0.75 
and T= 12800 K, however we consider the higer temperature as 
more probable. The IR SED can not be fitted with one BB. We 
derive two dust components, with 900 K and 420 K. 
Var A. In the period 1997 - 2004 the star did not change its bright- 
ness notably ('Humphreys et al. "2006). From our spectrum we esti- 
mate a possible temperature range as 7000 < T < 15000 K. Ac- 
tually there are only H I emissions and absorption blends of Fe II 
in the spectrum. A nearby nebula located in 9^^-12^^ to NE di- 
rection, gives Av = 0.9 — 1.1. This is the on ly star in our list 
whose luminosity was determined quite firmly jHumphrevs et al.l 



Figure 5. Spectral energy distributions of the studied stars in M33. 
There are observed data (squares), data corrected for interstellar absorp- 
tion (crosses), the best fit BB models of hot and dust components (dashed 
lines) and sum of the spectral components (solid lines). 

1 19871 : iHumphrevs eTaD l2006h. Ig L/Lq ^ 6.0 (with the distance 
to M33 adopted by us). These authors argue that the star has a 
notable amount of neutral extinction. To be consistent we use the 
same method as we do for other stars, but we consider this analy- 
sis as formal. We find that the luminosity IgL/LQ =6.0 may be 
reached at Ay = 2.5 (T=10000K). One needs better spectrum of 
VarA to find correct extinction and the star temperature. 
Var B. The star shows strong variability. After its max imum of V ^ 
15 mag in 1992 it weaken e d to 17.7 niag in 2 005 dSzeif ert et al.l 
| l99d lMassevetal.1 12OO6I : IViotti et al.l I2OO6I: IShporer & Mazeh| 
2W6). We used only UBV data from Viotti et al. (2006) and Spitzer 
data as about simultaneous. A surrounding nebula gives us the ex- 
tinction limits Av = 1.4 — 2.8. We estimate the star's temperature 
range as 15000 - 25000 K (strong He I, [Fell], [Felll] lines, and 
no Bowen blend or Hell). Fitting the SED, we found the maximal 
value of Av we may adopt to be 1.0 mag (corresponding to a tem- 
perature of 25000 K). Otherwise the te mperature is too high Note 
that in another brightness state of the star lSzeifert et in (ll996h used 
Av = 0.7. They estimated extinction as foreground galactic plus 
1/2 of the value estimated from HI measurements. In spite of the 
different methods the star luminosity they determined is the same 
as that we found. 

Var C. In 1985 the star showed a maximum brightness of B=15.4, 
and in 1997-2001 a minimum of V^17.3 with a variability of less 
than 0.2 mag, although it beg un to^brighten again until 2004 on 
^ 1 mag ( Szeif ert et aLlll996l : IShporer & Mazehll2006h . We used 
our main sources of the data. A nebula measured to SE direction 
from the star to 11^^ gives in different parts Av from 0.32 to 0.73. 
The star temperature estimate is 10000 - 15000 K. There are Mg II 
lines, bright hydrogen emissions, and no He I lines. Fitting of the 
SED gives best values of T= 14000 - 16000 K and Av = 0.4. At 
Av > 0.45 the temper ature becomes higher than we may accept 
from our spectrum. Note I Szeifert et al.l ( Il996h used Av — 0.8, and 
their luminosity derived is 10% higher. 
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Figure 6. Temperature - luminosity diagram for seven studied M 33 stars. 
In N 93351 and Var83 we used two most probable values of extinction, 
inner area shows the acceptable L, T values. In Var A our analysis is rather 
formal (see text). In other stars we used one value of interstellar extinction 
(Table 2 ). The segments indicate a range of acceptable temperatures. Stellar 
evolution tracks ( Claret 2004) for Z=0.007 stars are plotted. 



5 CONCLUDING REMARKS 

Fig.[6] shows the temperature - luminosity di agram for th e seven 
studied stars in M33. Stellar evolution tracks (IClaretll2004l) for 40 
-125 Mq stars are plotted, which were abridged in the latest evolu- 
tion stages so as not to crowd the figure. All these stars have masses 
of 40 — 130 Mq, and are located in the region of evolved LBV 
stars. The new star N 93351 is among the brightest evolved stars in 
M33. 

One finds from Table 2 that the dust reradiates on average 
8 % of the energy budget of the stars. The total scatter of stellar 
component luminosities is 1.0 dex, while that of the dust luminosi- 
ties is 2.0 dex. Any possible uncertainties in interstellar extinction 
are not as important in IR radiation as they are in the optical range, 
but we observe much bigger scatter in the dust component lumi- 
nosities than that in stellar components. This may mean that the 
dust phenomenon in the most massive evolved stars is quite tempo- 
rary, and that the dust components must be variable on scales much 
less than those of core evolution. One may conclude that the dust 
activity of LBV and related stars follows the strong eruptions (at 
time scale of tens, or p ossibly, hundreds of yea rs), as was detected 
in Var A in 1950-1953 (iHumphrevs et al.lll987h . 

The new star N 93351 has properties intermediate between 
the LBVs and Var A , in the same time it has a hot photosphere, 
LBV-like luminosity and an extensive curcumstellar material. This 
makes N 93351 very important for detail studies in future. 



Var 83. The star did not sho w strong variability in 1997 - 2004 
jMocheiska et al.l [200 ll : iMa ssev et al. 2006); its variations were 
less than 0.5 mag. Its nebula indicates Ay ~ 1.35. In our spec- 
trum there are very strong He I, Fe II, [Fe II], [Fe III] and Si II emis- 
sions (no Bowen blend), which implies a temperature range 15000 
- 25000 K. Our SED solutions give T> 30000 K at Ay ~ 1.35, 
which is in contradiction with the spectrum. Hence we use two val- 
ues Av = 0.7 (Fig.[5]) and 1.0 giving the m ost close temperatures 
to those estimated from the spectrum. Note ISzeifertetal1(ll996h 
used Av = 1.05, they obtained the star luminosity inside the range 
of our limits Table 2 . 

Var 2. The star has not shown strong variabilit y since 1935 (< 0.3 
mag), now its brig htness is about V ^ 18.2 (fszeifert et al.lll99d : 
iMassev et aLl2006h . There are no 2MASS and Spitzer data (the star 
is fain t). We show the K band point in Fig.|5]from Szeife rt et al.l 
( 1996) as an upper limit. We found Av ~ 0.85 in Var 2's nebula. 
There are He I, He II and Bowen blend lines in Var 2 spectrum. The 
temperature is estimated as 35000 - 50000 K. The best fit of the 
SED gives estimates Av ~ 0.8 and 40000 K. 
V 532. The star shows strong historical variab ility, however in the 
period 1997 - 2005 it was relatively calm (iKurtev et all 1200 ll : 
IZharova & Sholukhoval [2004 IViotti et aLll200l and our own un- 
published data), from B ^ 17.2 in 1997 to 17.4 in 2000, drop- 
ping to 17.7 to October 2001, and returning to 17.4 in 2004 - 
beginning of 2005. A close nebula surrounding V532 indicates 
Av 1.1. There are Hell, Bowen blend and strong He I emis- 
sions in our spectra, which implies that the temperature has to be 
30000 - 40000 K. The best fit of the SED we find with T = 35000 K 
(Av = 0.8). 
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